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ABSTRACT 


The  purpose  of  this  study  is  to  develop  a  computer  program  that  will  predict  the 
damage  progression  in  composite  plates  subjected  to  bending  loads.  Kwon’s 
micromechanical  model  is  used  to  compute  the  smeared  effective  moduli  from  the  material 
properties  of  fiber  and  matrix  as  well  as  to  determine  stresses  at  the  constituent  level. 
Failure  criteria  based  on  micro-stresses  are  then  applied  to  determine  the  extent  and  type  of 
damage  that  occur  in  the  composite  under  various  loading  conditions.  The  progression  of 
damage  throughout  the  composite  until  complete  failure  of  the  composite  can  then  be 
simulated  using  the  current  computer  program.  The  numerical  prediction  for  a  laminated 
composite  plate  containing  a  hole  and  subjected  to  a  bending  load  agrees  well  with  the 
experimental  data. 
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I.  INTRODUCTION 


Composite  materials  are  becoming  more  competitive  with  metallic  materials  in 
structural  applications.  Fiber  reinforced  composites  are  the  materials  of  choice  when  design 
factors  dictate  the  need  for  high  strength  to  weight  ratio,  high  stiffness  to  weight  ratio  and 
resistance  to  harsh  environmental  conditions.  Aircraft  and  aerospace  industries  have  in 
particular  taken  advantage  of  fiber  reinforced  composite  to  decrease  weight  while 
maintaining  strength.  Naval  applications  of  composites  are  increasingly  more  common  as 
they  have  become  more  of  an  economically  viable  option.  Most  composite  structures  are 
usually  subjected  to  bending  loads.  It  is  therefore  desirable  to  predict  the  strengths  and 
failure  modes  of  composite  plates  subjected  to  bending  loads  such  that  full  advantage  of  the 
plates’  material  properties  may  be  used. 

Predicting  failure  mechanisms  of  composite  structures  is  more  difficult  than  for 
metallic  structures.  In  addition,  most  studies  dealing  with  composite  failure  have  been  with 
inplane  loading  of  composite  plates.  It  is  the  purpose  of  this  study  to  develope  an  analytical 
model  to  determine  the  failure  modes  and  strength  of  laminated  fibrous  composite  plates 
subjected  to  bending  loads. 

Chang  and  Chang  studied  failure  in  a  laminated  composite  plate  containing  a  hole 
at  its  center.  The  plate  was  subjected  to  tensile  loading.  A  finite  element  program  was 
created  to  model  the  experimental  results.  Stress  analysis  and  failure  analysis  were 
addressed  separately  in  the  program.  The  computer  program  modeled  stress  based  on  the 
classical  lamination  theory  of  composites.  The  plate  load  was  increased  incrementally  so 
that  a  progressive  damage  study  of  the  plate  could  be  performed.  As  elements  failed,  their 
material  properties  were  degraded  by  the  computer  program  so  that  failed  elements  would 
no  longer  carry  any  of  the  applied  load.  The  elements  surrounding  failed  elements  then 
assumed  the  load  previously  supported  by  the  failed  elements  and  experienced  increased 
stress  due  to  those  loads.  (Chang,  1987) 

The  failure  analysis  was  based  on  work  done  by  Yamada  and  Sun  (Yamada,  1978). 
Matrix  cracking,  fiber  and  matrix  shearing,  and  fiber  breakage  were  the  identified  modes 
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of  failure  for  the  composite  constituents.  When  stresses  found  in  the  stress  analysis  exceed 
the  failure  criteria,  that  element  would  fail  and  damage  progression  in  the  composite  would 
be  computed.  The  computer  model  developed  by  Chang  and  Chang  was  determined  to 
provide  an  accurate  predication  of  composite  damage  propagation.  (Chang,  1987) 

Owen  and  Li  stated  that  for  composite  plates  a  macromechanical  approach 
involving  the  use  of  the  classical  laminated  plate  theory  was  not  sufficiently  accurate  for 
stress  calculations,  yet  a  micro  approach  would  be  too  numerically  cumbersome  to  provide 
any  results.  Laminated  composite  plates  have  been  shown  to  have  the  greater  effect  of 
tran.sver.se  shear  deformation  as  the  plate  thickness  increases  relative  to  the  plate  length  and 
as  such  the  classical  laminated  plate  theory  is  not  sufficiently  accurate  to  model  orthotropic 
and  anisotropic  laminates.  Owen  and  Li  used  a  model  based  on  inplane  displacements 
across  the  plate  thickness,  in  which  inplane  displacements  are  pieces-wise  linear  and  lateral 
displacements  constant  across  the  thickness.  This  allowed  for  the  three  dimensional  model 
to  be  reduced  to  a  two  dimensional  plate  bending  model.  (Owen,  1987) 

Kwon’s  use  of  micromechanical  analysis  provides  the  means  to  determine  damage 
progression  at  the  constituent  level  for  laminated  fibrous  composites.  The  micromechanical 
model  is  used  to  determine  smeared  composite  properties  from  the  fiber  and  matrix 
material  properties.  Macro  analysis  is  then  used  to  determine  displacements  and  strains  of 
composite  using  the  smeared  composite  properties  in  the  finite  element  analysis. 
Microstresses  are  determined  for  each  composite  constituent  in  the  element  by  using  the 
micromechanical  model  in  conduction  with  the  composite  strains.  Failure  analysis  can  then 
be  applied  at  the  fiber  and  matrix  level.  Damage  progression  of  the  composite  may  be 
predicted  based  on  the  strengths  and  failure  modes  of  the  fibers  and  matrix  that  comprise 
the  laminated  composite.  (Kwon,  1993) 
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II.  MICROMECHANICAL  MODEL 


A.  OVERVIEW 

The  investigation  of  damage  progression  in  a  laminated  fibrous  composite  plate  is 
undertaken  using  a  micromechanical  model.  Kwon’s  micromechanical  model  permits  the 
calculation  of  microstresses  for  a  composite  using  the  material  properties  of  the  fiber  and 
matrix  that  make  up  that  composite  (Kwon,  1993).  The  microstresses  are  evaluated  with  a 
set  of  failure  criteria,  and  thus  a  study  of  damage  throughout  the  composite  with  details  of 
fiber  and  matrix  damages  is  possible.  The  mircromechanical  model  also  computes  the 
smeared  composite  material  properties  from  the  fiber  and  matrix  material  properties. 

B.  MICROMECHANICAL  CELL 

The  micromechanical  unit  cell  is  made  of  a  fiber  and  surrounding  matrix  material. 
Using  the  symmetry  of  the  cell,  only  a  quarter  of  the  actual  cell  needs  to  be  modelled  as 
shown  in  Figure  1 .  The  fiber  direction  is  taken  as  the  direction  along  the  1  axis.  The  cell  is 


Figure  1.  Kwon’s  Unit  Cell  and  Subcells 

divided  into  four  subcells.  Each  subcell  size  depends  on  the  volume  fraction  of  the  fiber,  ¥. 
(Kwon,  1993) 
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The  stresses,  ay,  and  strains,  Ey,  of  the  composite,  are  the  volume  average  of  the 
subcell  stresses,  a“y,  and  subcell  strains,  e“y,  as  given  below  (Kwon,  1993); 


/a" 


+  v/  1 


J7) 


a  ij  + 


1  -  a/7 ]a./  +  f  1  -  77l  a,.,. 


Eq  1 


E..  = 
u 


/e"  +  /(  1  -  77  ]e*,y  +  /[  1  -  77  )e  /  +  [  1  -  77 
where  a  is  the  fiber  subcell,  b,  c,  and  d  are  the  matrix  subcells. 


Eq2 


Stress  continuity  at  subcell  interfaces  is  satisfied  by  the  following  relationships 
(Kwon,  1993): 


a  c  da  c  b  d 

a  22  =  Ct  22,  O  22  =  O  22,  a  33  =  O  33,  033  =  033 

b  c  da  c  h  d 

<7  12  =  ct  12,  a  12  =  o  12,  a  13  =  a  n,  0,3  =  aj3 


—(‘bed 
Cf  23  =  O  23  =  Cy  22  =  O  23 


Eq3 

Eq4 


The  strains  of  the  subcells  are  assumed  to  meet  the  following  requirements  (Kwon, 


1993): 


a  b 

e  11  =  e  11 

a  h 

e  22  +  e  22 

a  c 

e  33  +  e  33 


c  d 

£  11  =  e  11 

c  d 

e  22  +  e  22 

b  d 

e  33  +  e  33 


Eq5 


abed 

e  12  +  e  12  =  e  12  +  e  12 

a  c  b  d 

e  13  +  e  13=  e  13  +  e  13 


Eq6 


The  constitutive  equation  for  each  subcell  based  on  Hooke’s  law  is  given  by; 

Eq7 


«  _u  u 

O  ij  =  E  ijkiE  kl 


where  i,  j  =  1 , 2,  3  and  a  =  a,  b,  c,  d. 

The  stre.ss  for  each  subcell  expressed  in  Equation  7  is  substituted  into  Equation  3 
and  Equation  4  to  have  an  expression  for  stress  continuity  represented  in  terms  of  subcell 
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strains.  The  resultant  equations  along  with  Equation  2,  Equation  5  and  Equation  6  are 
solved  simultaneously  to  result  in  subcells  strains  given  by  average  composite  strains. 
Subcell  strains  expressed  in  terms  of  the  average  composite  strains  are  substituted  into 
Equation  7  and  the  resultant  subcell  stresses  expressed  in  terms  of  the  composite  strains  are 
plugged  into  Equation  1  to  result  in  the  relationship  between  the  composite  stresses  and 
composite  strains.  This  relationship  gives  the  smeared  composite  material  properties 
expressed  in  terms  of  the  fiber  and  matrix  material  properties.  (Kwon,  1993) 

As  a  result,  the  micromechanical  model  computes  the  composite  material  properties 
using  the  constitutive  fiber  and  matrix  material  properties.  In  addition,  the  model  gives  the 
relationship  between  microstresses  (and  strains)  and  composite  stresses  (and  strains).  That 
is,  microstrains  are  found  from  composite  strains.  The  microstrains  are  then  used  with 
Hooke’s  law  for  each  constituent  to  determine  the  microstresses,  from  which  the  composite 
stresses  are  determined. 
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III.  FINITE  ELEMENT  MODEL 


A.  PLATES  IN  BENDING 

One  objective  of  this  study  is  to  develope  a  computer  program  that  will  model  the 
failure  experienced  by  a  thin  composite  plate  due  to  a  bending  load.  The  finite  element 
method  along  with  Kwon’s  micromechanical  model  is  used  for  the  present  study.  The  plate 
bending  element  used  in  this  study  is  described  below. 

A  plate  bending  element  is  shown  in  Figure  2  where  x,  y,  and  z  describe  the  global 


coordinates  of  points  of  the  plate  and  u,  v,  and  w  are  the  displacements  of  those  points,  t  is 
the  plate  thickness.  The  xy  plane  is  parallel  to  the  midsurface  plane  prior  to  deflection. 
The  displacement  of  any  point  in  the  plate  can  be  expressed  as  (Ugural,  1987). 

{d}  ^  =  {u(x,y,z),vi  (x,  y,z),w  (x,  y) ) }  ^  Eq  8 

That  is,  the  inplane  displacements  vary  through  the  plate  thickness  while  the  transverse 
displacement  remains  constant  through  the  plate  thickness.  This  displacement  within  the 
element  may  be  written  in  terms  of  shape  functions  and  nodal  displacement. 

n 

=  Eq9 

/ 

Equation  9  can  be  rewritten  in  a  matrix  form  as  given  below  (Owen,  1987), 

{d}^,  =  [/V]  (61^  EqlO 
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where  {6}^  is  the  nodal  displacement  vector  for  a  plate  element.  {5};  is  the  nodal 
displacements  for  the  i'*’  node  of  an  element. 


15},= 


Eq  11 


where  superscripts  b  and  t  denote  the  inplane  displacements  at  the  bottom  and  top  sides  of 
the  plate  element,  respectively.  A  four  node  quadrilateral  element  is  used  for  the  present 
•Study  for  interpolation  of  displacements  in  the  x-y  plane  as  shown  in  Figure  3; 


4 


2 


i  s 


1 

r 

j)  3 

^  7 

(0,0) 

(b - ( 

)  1 

Figure  3.  Quadrilateral  Element 


The  shape  function  for  the  quadrilateral  element  is  given  in  terms  of  the  natural  coordinate 
.system: 


/V,  =  ^(l+/T,)  (1+.SW.)  Eql2 

where  i=  1,2, 3, 4,  and  q  and  sj  are  the  local  coordinates  of  the  corners  of  the  quadrilateral 
element  whose  values  are  either  -1  or  +1. 
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Linear  shape  functions  are  given  below  in  terms  of  the  third  natural  coordinate  axis, 
which  is  mapped  into  the  plate  thickness  direction  in  the  physical  coordinate  system: 


Eq  13 


Combining  the  previous  bilinear  shape  functions.  Equation  12  with  the  linear  shape 
functions.  Equation  13,  yields  (Kwon,  1988): 


Nij-x  =N.{r,s)H^{q)  , 

N2j  =  Nj{r,s)H2{q) 


1,2,  3,4 


Eq  14 


Three  cartesian  displacements  within  a  plate  element  are  then  written  as: 

=  N lU  \  +  N 2U  \  +  N^U  2  +  ^4^  2  +  3  +  N(,U  3  +  N-jU  4  +  N^U  4 

=  N^V^I  +  N2v‘i  +  Nt,V^2  +  ^4^2  +  +  A'7V^4  +  A^gV*4 

^2^1  ^4*^4 

where  u  and  v  are  inplane  displacements  and  w  is  the  transverse  displacement.  Subscript 
and  superscript  for  the  nodal  displacements  indicate  the  node  number,  and  top  or  bottom 
plane  displacement  of  the  plate,  respectively.  Since  there  is  no  rotational  degree  of  freedom 
in  the  present  plate  bending  element,  the  inplane  displacements  have  two  components  at 
one  nodal  point:  i.e.  top  plane  displacement  and  bottom  plane  displacement. 

When  a  bending  load  is  applied  to  a  plate  to  cause  it  to  deflect  from  its  original 
configuration,  deflection,  stresses  and  strains  associated  with  the  deformation  will  follow 
several  basic  assumptions.  These  assumptions  are  (Ugural,  1987): 

1)  The  normal  stress  to  the  midsurface,  a^,  is  negligible 

2)  Deflection  of  the  plate  is  small  in  comparison  with  the  thickness  of  the  plate 

3)  The  normal  strain  in  the  plate  thickness  direction,  may  be  neglected. 
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Strains  for  the  element  are  in  a  matrix  form: 


^xy 


lyz 


1 A 

d_y  dx 


0 

0 

0 


dz 

0 


A  A 

dz  dy 


u 

V 


w 


e 


Eq  16 


Equation  10  is  used  to  represent  the  displacements  in  Equation  16  so  that  strains  can 
be  expressed  in  terms  of  the  shape  functions  and  displacements.  Elemental  strain  in  terms 
of  shape  functions  and  displacements  are  shown  in  the  Appendix.  The  matrix  of  shape 
function  derivatives  is  referred  to  as  [B].  Strain  equations  may  be  expressed  as: 

=  Ifi]  {5}^  Eq  17 

Stresses  of  a  plate  in  accordance  with  the  assumptions  listed  for  this  study  are: 


{(?} 


Eq  18 


Stres.ses  are  obtained  from  the  strains  using  Hooke’s  law.  For  an  isotropic 
homogenous  plate,  stres.ses  are  (Ugural,  19X7): 


10 


1  2  ,  2 
1  -V  1  -V 


2  ,  2 
-V  1  -V 


0  0 
0  0 
0  0 


E 

2(1  +v) 


E 

2(l+v) 


E 

2(l+v) 


where  E  is  the  elastic  modulus  and  v  is  poisson’s  ratio. 

The  matrix  of  material  properties  may  be  written  in  a  more  condensed  form  and 
called  [D],  This  allows  Equation  19  to  be  written  as: 

{CJ},  =  [d]  {£},.  Eq20 

For  a  unidirectional  fibrous  composite,  the  [D]  matrix  is: 


Djj  £>12  0  0  0 

D12D22  0  0  0 

0  0  D33  0  0 

0  0  0  D44  0 

0000 


”  “  1 

^12^21 

EjjH2i 

12  1 

-  ^112^21 

^22 

22  1 

-  M^12M^21 

D33 

=  Cj2 

=  Gi3 

^55 

“  ^23 
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in  which  Ej  j  is  Young’s  modulus  in  the  direction  of  the  1  axis  as  shown  in  Figure  1 ,  E22  is 
modulus  along  the  2  axis,  H12  is  poisson’s  ratio  on  the  1  face  in  the  2  direction,  |i2i  i-s 
possion’s  ratio  on  the  2  face  in  the  I  direction,  and  Gjj  is  the  shear  modulus  associated  with 
i  and  j  axes. 

Equation  20  shows  that  stresses  are  obtained  from  strains  and  those  strains  are 
determined  using  displacements  in  Equation  17.  A  method  now  must  be  developed  to 
determine  displacements  such  that  stress  can  be  obtained.  Energy  methods  can  be  used  to 
find  displacements.  The  Rayleigh-Ritz  method  is  particular  suited  to  determining  unknown 
displacements.  The  procedure  of  using  the  Rayleigh-Ritz  method  is  (Ugural,  1987): 

1 )  Assume  a  displacemnet  function  in  terms  of  unknown  nodal  parameters 

2)  Compute  the  potential  energy,  O,  in  terms  of  the  unknown  parameters  which 
describe  the  total  work,  W,  and  strain  energy,  U 

n  ^  U-W  Eq23 

3)  Apply  the  principle  of  minimum  potential  energy  and  solve  for  the  displace¬ 
ments. 

Strain  energy,  U,  may  be  expressed  in  terms  of  stress  and  stain  (Ugural,  1987). 


where  V  is  the  volume  of  the  body. 

The  work,  W,  done  by  external  forces,  F  (per  unit  volume),  and  surface  forces,  T 
(per  unit  area)  is  given  by  (Ugural,  1987): 

^  =  j  ( ^  dV  +  T^u  +  T^.v  +  w)  dS  Eq  25 

V'  A 

where  A  is  the  boundary  surface  of  the  body. 

Using  the  principle  of  minimum  potential  energy,  fl  is  minimized  such  that 
(Ugural,  1987): 

Sn  =  8  (a,e^  +  +  x„y,y  +  dV  E, 26 

-f \^F^ui-F^v  +  F^w)dV +\{T^u  +  T  v+T^w)ds'^  =  0 
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The  potential  energy  equation  can  be  written  using  the  finite  element  discretization 


n  =  {d]J{F}^)dV-Y,\ {d}J{T}^dS  Eq27 

where  summations  are  over  the  discretized  finite  elements,  and  {d}  is  the  displacement 
vector.  Stresses  in  Equation  27  can  be  written  in  terms  of  strains  and  [D]  from  Equation  20. 

n  =  <•=>.-  Eq28 

1  s 

Displacements  for  each  element  are  represented  by  the  unknown  nodal 
displacements  {5}„  and  shape  functions.  Strains  can  now  be  expressed  in  terms  of 
elemental  nodal  displacements  and  derivatives  of  the  shape  functions  (Owen,  1987).  Using 
Equation  17  for  the  strain  in  terms  of  the  [B]  matrix  and  nodal  displacements,  and  Equation 
10  for  the  displacement  function.  Equation  28  becomes: 

n=  [B]  {5}^-{5}^[/V]^{F)^]rfU  Eq29 

-2|  {b}^[Nf{T]^dS 

s. 

After  minimizing  the  total  potential  energy  with  respect  to  the  nodal  displacements, 
the  element  stiffness  matrix  is  defined  as  (Owen,  1987): 

[k],=  J  [Bf[D]  [B]dV  Eq30 

The  element  force  matrix  is  formed  by  combining  the  work  done  by  external  surface  forces 
and  body  forces: 

{<21,  =  J  [Nf{F}dV  +  j[Nf{T]dS:  Eq31 

Assembling  all  element  matrices  and  vectors,  the  system  equation  is: 

[AT]  {5}  =  {Q}  Eq32 

The  unknown  nodal  displacement  can  thus  be  determined. 
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Since  the  shape  of  the  elements  can  be  irregular  quadrilateral,  domain  and  boundary 
integrations  of  such  shapes  are  solved  using  the  natural  coordinate  system  (r,  s,  q) 
(Segerlind,  1984).  This  integration  must  be  done  numerically  and  will  use  the  Gauss- 
Legendre  integration  technique.  Since  the  natural  coordinate  system  is  used  for  numerical 
integration  to  compute  the  integral,  the  Jacobian  matrix,  [J],  must  be  used  for  mapping 
between  the  physical  and  natural  coordinate  systems,  that  is. 


[J] 


dx  By  Bz 
Br  Br  Br 

Bx  By  Bz 
Bs  Bs  Bs 

Bx  By  Bz 
Bq  Bq  Bq 


Eq33 


The  stiffness  matrix  can  now  be  represented  in  terms  of  the  natural  coordinate 
system  (Segerlind,  1984). 


1  I  1 

[^],  =  J  [5]^[D]  [B]dV  =  J  J  J  [fi]^[Z)]  [B]\J\drdsdq  Eq34 


The  integrals  for  the  natural  coordinate  system  are  defined  on  the  interval  +1  to  -1. 
Using  numerical  integration  with  n,  m,  1  integration  points  along  each  direction,  the 
.stiffne.ss  matrix  becomes: 


[A],  =  \[B\‘  [D\  \B\dV 

V 


m  n  l  , 

I  Z  Z  [B||y| 


i  =  \j=  \  k=  \ 


where  indicates  the  function  is  evaluated  at  (q,  Sj,  q^).  That  is,  [fi],  [D],  and  L/l  are 

evaluated  at  their  sampling  points  and  W  is  the  weight  coefficients  for  Gauss  quadrature 
(Segerlind,  1984). 


With  the  calculation  of  a  stiffness  matrix  and  a  force  vector,  displacements  for  the 
nodal  points  of  the  plate  can  be  determined.  These  displacements  are  used  with  Equation 
17  to  solve  for  the  strain.  Once  the  strains  are  determined,  microstrains,  microstresses  and 
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macrostresses  are  determined  using  the  micromechanical  model  described  in  the  previous 
chapter. 
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IV.  MICRO-MACRO  ANALYSIS 


A.  MICRO-MACRO  ANALYSIS 


To  find  the  damage  progress  in  a  plate  due  to  an  external  loading,  a  micro¬ 
macromechanical  approach  is  used.  Kwon’s  micromechanical  model  is  the  basis  for  the 
procedure.  Composite  material  properties  are  determined  from  the  fiber  and  matrix 
material  properties  using  the  micromechanical  model.  These  composite  properties  are  then 
used  to  determine  nodal  displacements  from  the  finite  element  analysis  (macro  analysis). 
Equation  17  shows  that  macro-strains  are  obtained  from  displacements.  The 
micromechanical  model  is  again  used  with  the  macro-strains  to  derive  the  micro-strains  of 
each  element.  Micro-stresses  are  calculated  using  micro-strains.  Failure  criteria  applied  to 
measure  damage  progression  in  the  composite. 

[B]  can  be  decomposed  into  two  parts;  [B]f,for  the  bending  strains  and  [B]^  for  the 
shear  stains  (Kwon,  1994). 


du 

dx 

dv 

dy 

I 

du  dv 
_dy  ^  dx_ 


Eq35 


shear  (/) 


du  dw 
dz  dx 
5v  bw 
5z  ^  5y 


Eq36 


[B]i3  is  evaluated  at  four  integration  points  in  the  rs  plane  while  [B]^  is  evaluated  at  one 
integration  point  to  prevent  shear  locking.  As  a  result,  the  bending  strains  are  averaged  to 
produce  one  value  for  each  element  in  the  same  manner  as  the  transverse  shear  strains. 
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B. 


FAILURE  CRITERIA 


In  order  to  study  the  damage  progression  in  a  composite  subjected  to  a  bending 
load,  the  following  failure  criteria  are  used  (Berner,  1993); 


o 


> 


xj 


Eq37 


O22 


a 


m\ 


22 


V  X', ) 


O12 


V  s'"  j 


m\2 

!i2_ 

V  s'"  j 


>  1  if  o''‘22  >  0 

>  1  if  a"‘22  <  0 


Eq38 


in  which  a*^|  j  is  the  longitudinal  fiber  stress  and  is  the  ultimate  strength  of  the  fiber  in 
the  longitudinal  direction.  Subscripts  t  and  c  denote  tensile  and  compressive  strengths. 
a'^22  is  the  transverse  matrix  stress,  X'"  is  the  ultimate  strength  of  the  matrix,  a"’ 12  is  the 
matrix  shear  stress,  and  is  the  ultimate  shear  strength  of  the  matrix. 

The  failure  criteria  are  used  in  determining  the  failure  of  the  composite  in  four 
different  modes;  fiber  failure  in  compression  or  tension  and  matrix  failure  in  compression 
or  tension.  If  any  of  these  failure  modes  occurs,  then  the  material  properties  for  the  fiber 
and/or  matrix  related  to  the  failure  are  degraded  at  the  failed  location.  This  degradation  of 
material  properties  in  the  computer  modeling  is  accomplished  by  multiplying  those 
material  properties  by  0.005,  effectively  reducing  the  properties  to  small  values  but  non¬ 
zero  values  allowing  the  material  property  matrices  to  be  non-singular.  For  subsequent 
calculations,  those  failed  fiber  and/or  matrix  will  not  take  any  of  the  applied  load.  (Since 
the  failed  material  properties  are  not  exactly  zeros,  the  failed  location  may  still  carry  very 
small  insignificant  loads.) 

The  load  is  applied  incrementally.  An  initial  load  is  chosen  and  applied  to  the 
composite.  The  computer  program  checks  whether  computed  stresses  satisfy  the  failure 
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criteria.  If  there  is  no  failure  of  either  the  fiber  or  matrix  in  the  composite,  the  program 
increases  the  load  incrementally  as  directed  by  the  input  file  until  some  type  of  failure 
occurs  or  the  maximum  prescribe  load  is  applied. 

If  some  type  of  failure  has  occurred,  then  that  fiber  and/or  matrix  material  is 
degraded  as  described  previously.  Another  iteration  will  now  be  carried  out  using  this 
“damaged”  element(s).  Each  loading  condition  is  allow  to  run  through  only  five  iterations. 
This  is  done  to  more  closely  model  actual  loading  of  the  plate  by  a  testing  machine  that 
continually  increases  the  load  over  a  finite  period  of  time.  The  failure  of  the  composite 
itself  was  based  on  the  degree  and  direction  of  the  failure  for  the  fibers  and  matrix. 

The  failure  is  checked  at  the  numerical  integration  points  of  each  plate  bending 
element.  The  elemental  stiffness  matrix  is  the  sum  of  stiffness  matrices  for  bending  and 
shear  (Kwon,  1994). 

+  Eq39 

For  the  bending  stiffness,  2-point  integration  rule  is  used  in  every  direction.  However,  1- 
point  integration  rule  is  used  for  the  shear  stiffness  on  the  inplane  directions  and  2-point 
rules  on  the  thickness  direction.  As  a  result,  the  stresses  due  to  bending  are  computed  at 
four  points  in  the  plane  of  a  plate  element  while  the  transverse  shear  stresses  at  one  point. 
The  bending  stresses  are  averaged  for  four  points  on  both  tension  and  compression  sides  of 
a  plate  element.  Finally  the  averaged  bending  stresses  and  the  shear  stresses  are  used  for 
the  failure  criteria  to  represent  failure  of  the  element  either  on  the  tension  or  compressive 
side. 
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V.  RESULTS 


A.  PLATE  WITH  A  HOLE  SUBJECTED  TO  UNIFORM  BENDING 

A  numerical  prediction  was  conducted  using  the  present  micro-macro  damage 
model  for  the  laminated  composite  used  by  Yang  in  his  experimental  testing.  The  laminated 
composite  plate  was  a  carbon/epoxy  composite  (IM6/3501-6)  manufactured  by  Hercules 

Materials  Company  of  Magna,  Utah.  The  fibers  were  in  aligned  in  the  0®  direction.  Yang 
used  a  12  inch  by  1  inch  plate  containing  a  0.250  inch  diameter  hole  at  the  center  and 
subjected  to  a  four  point  bending  test  to  create  an  uniform  bending  moment  along  the  edge 
of  the  plate.  (Yang,  1994) 

The  material  properties  of  the  composite  used  in  the  bending  study  were  determined 
from  experiments.  (See  Table  1) 


m 

Elastic  Modulus(Msi) 

2.30 

0.200 

Poisson’s  Ratio 

0J70 

0.360 

Fiber  volume  fraction 

0.250 

N/A 

Table  1.  Material  Properties  of  Composite 
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Yang  showed  that  for  the  one-quarter  symmetric  sample  of  the  plate,  as  in  Figure  4, 
the  failure  loading  corresponded  to  230  lb-in  applied  moment. 


Finite  Element  Grid  for  Symmetric  Plate  with  Hole 


Figure  4.  Symmetric  Plate  with  Hole 


Failure  of  each  element  of  the  plate  is  based  on  the  criteria  of  Equation  37  and 
Equation  38.  Total  failure  of  the  composite  is  determined  from  when  there  is  nearly 
complete  failure  of  the  fibers  along  the  y-axis.  Matrix  failure  occurs  in  different  forms. 
Matrix  may  fail  along  the  y-axis  creating  a  matrix  cracking  and/or  it  can  fail  along  the  x- 
axis  to  cause  fiber  splitting.  The  underlying  definition  of  composite  failure  in  this  case  is 
the  inability  of  the  composite  to  carry  an  applied  load.  The  damage  progression  predicted 
using  the  present  micro-macro  model  is  plotted  in  Figure  5  through  Figure  8.  As  seen  in 
Figure  5  and  Figure  6,  matrix  cracking  occurs  along  the  fiber  direction  causing  fiber 
.splitting.  As  the  load  increa.ses,  fiber  breakage  occurs  along  the  minimum  section  of  the 
plate  until  the  final  fracture  at  200  lb-in  the  failure  modes  agree  with  the  experimental 
ob.servation  and  the  final  failure  load  is  15%  le.ss  than  the  experimental  failure  load. 
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Fiber  Failure  Compression  Strength  Equals  Tensile  Strength  with  94.00  lb-in  Applied  Moment 


r- 

00 

d 


<N 

VO 

d 


(D 

o  o 


’Tt 

d 


S^ljOU] 


Figure  5.  Fiber  Failure  Compression  Strength  Equal  to  Tensile 
Strength  with  94.00  lb-in  Applied  Moment 
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Fiber  Failure  Compression  Strength  Equals  Tensile  Strength  with  165  lb-in  Applied  Moment 
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Figure  6.  Fiber  Failure  Compression  Strength  Equal  to  Tensile 
Strength  with  165  lb-in  Applied  Moment 
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Fiber  Failure  Compression  Strength  Equals  Tensile  Strength  with  188.00  b-in  Applied  Moment 
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Figure  7.  Fiber  Failure  Compression  Strength  Equal  to  Tensile 
Strength  with  188.00  lb-in  Applied  Moment 
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Fiber  Failure  Compression  Strength  Equals  Tensile  Strength  with  200  lb-in  Applied  Moment 
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Figure  8.  Fiber  Failure  Compression  Strength  Equal  to  Tensile 
Strength  with  200  lb-in  Applied  Moment 
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A  parametric  study  was  performed  by  varying  the  fiber  and  matrix  compressive 
strengths.  The  fiber’s  compressive  failure  strength  was  varied  while  all  other  composite 
material  properties  were  held  constant  as  before.  The  purpose  of  this  study  was  to 
determine  how  the  failure  load  of  the  composite  plate  varied  in  accordance  with  the  change 
of  the  compressive  strength  of  the  fiber.  The  matrix  compressive  strength  was  then  varied 
while  all  other  material  properties  of  the  composite  were  held  constant  to  analyze  its  effect 
on  the  failure  strength  of  the  composite.  Figure  9  and  Figure  10  show  the  reduced 
composite  strength  for  reduced  fiber  strength  and  increased  matrix  strength. 

In  cases  where  the  fiber  compressive  failure  strength  is  varied,  Figure  9  shows  that 
as  the  compressive  strength  of  the  fiber  was  decreased,  so  did  the  load  required  to  cause 
failure  of  the  composite.  The  reduction  in  the  required  load  is  approximately  equal  to  the 
percent  reduction  of  the  compressive  failure  strength  of  the  fiber.  A  straight  line  is  drawn 
in  Figure  9  for  the  cases  where  both  tensile  and  compressive  strengths  of  the  fiber  are 
reduced  uniformly.  The  composite  failure  load  for  the  reduced  fiber’s  compressive  strength 
only  is  larger  than  that  for  both  reduced  tensile  and  compressive  strengths  of  the  fiber,  as 
expected.  However,  there  is  no  more  than  a  5%  difference  between  the  two  composite 
failure  loads.  As  the  fiber  and  compressive  strength  is  reduced  more,  the  composite  failure 
strength  deviates  further  from  the  straight  line  because  the  matrix  strength  plays  a  more 
important  role  in  the  composite  failure.  Although  other  material  properties  may  influence 
the  finial  failure  load,  the  fiber  strength  is  the  dominant  material  property  in  determining 
failure  load. 

Figure  1 1  through  Figure  14  show  the  damage  progression  in  a  composite  for 
various  fiber  compressive  strengths.  The  results  from  these  figures  are  used  to  generate 
Figure  9.  As  shown  by  Figure  1 1  through  Figure  1 4,  the  applied  moment  required  for 
composite  failure  increases  as  the  fiber  compressive  strength  increases.  Composite  failure 
for  these  cases  initiates  in  the  fiber  splitting  mode  in  the  longitudinal  direction  and  fails  in 
the  transverse  fiber  breakage. 

Figure  10  is  the  graph  of  resultant  failure  loads  due  to  changing  the  matrix’s 
compressive  strength  while  keeping  all  other  composite  properties  the  same  as  listed  in 
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Table  1 .  Changing  the  matrix  strength  had  no  effect  on  the  fiber  failure  load  as  indicated  in 
Figure  14  through  Figure  23  which  show  the  damage  progression  in  a  composite  with 
various  matrix  compressive  strengths.  Matrix  failure  occurs  at  larger  bending  loads  for 
increased  matrix  strength  as  expected.  However,  the  composite  failure  load  changes  little 
for  increased  matrix  strength  because  the  fiber  is  the  dominant  material  for  ultimate  failure 
of  the  composite  structure. 
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Moment  Load  (lb-in) 


Percent  Compressive  Fiber  Strength  of  Ultimate  Fiber  Tensile  Strength 


Figure  9.  Failure  Load  with  Changing  Compressive  Fiber  Strength 
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Figure  10.  Failure  Load  with  Changing  Compressive  Matrix  Strength 
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Figure  1 1.  Fiber  Failure  Compression  Strength  Reduced  to  20%  of 
Tensile  Strength  with  47.00  lb-in  Applied  Moment 
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Fiber  Failure  Compression  Strength  Reduced  to  20%  of  Tensile  Strength  with  58.75  lb-in  Applied  Moment 
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Figure  12.  Fiber  Failure  Compression  Strength  Reduced  to  20% 
of  Tensile  Strength  with  58.75  lb-in  Applied  Moment 
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Fiber  Failure  Compression  Strength  Reduced  to  60%  of  Tensile  Strength  with  1 1 7.750  Ib-in  Applied  Moment 
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Figure  13.  Fiber  Failure  Compression  Strength  Reduced  to  60% 
of  Tensile  Strength  with  117.50  Ib-in  Applied  Moment 
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Fiber  Failure  Compression  Strength  Reduced  to  60%  of  Tensile  Strength  with  141.00  lb-in  Applied  Moment 
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Figure  14.  Fiber  Failure  Compression  Strength  Reduced  to  60%  of 
Tensile  Strength  with  141.00  lb-in  Applied  Moment 
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Matrix  Failure 
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Figure  15.  Matrix  Failure  Compression  Strength  Increased  to  200%  of 
Tensile  Strength  with  117.50  lb-in  Applied  Moment 
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Matrix  Failure  Compression  Strength  Increased  to  200%  of  Tensile  Strength  with  165  lb-in  Applied  Moment 


Figure  16.  Matrix  Failure  Compression  Strength  Increased  to  200% 
of  Tensile  Strength  with  165  lb-in  Applied  Moment 
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Matrix  Failure  Compression  Strength  Increased  to  200%  of  Tensile  Strength  with  188.00  lb-in  Applied  Moment 


Matrix  Failure  Compression  Strength  Increased  to  200%  of  Tensile  Strength  with  200  lb-in  Applied  Moment 
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Figure  18.  Matrix  Failure  Compression  Strength  Increased  to  200% 
of  Tensile  Strength  with  200  lb-in  Applied  Moment 
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Matrix  Failure  Compression  Strength  Increased  to  600%  of  Tensile  Strength  with  1 17.5  lb-in  Applied  Moment 
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Figure  19.  Matrix  Failure  Compression  Strength  Increased  to 
600%  of  Tensile  Strength  with  117.50  lb-in  Applied  Moment 
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Matrix  Failure  Compression  Strength  Increased  to  600%  of  Tensile  Strength  with  165  lb-in  Applied  Moment 
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Figure  20.  Matrix  Failure  Compression  Strength  Equal  to  600% 
of  Tensile  Strength  with  165  lb-in  Applied  Moment 
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Figure  21.  Matrix  Failure  Compression  Strength  Increased  to  600%  of 
Tensile  Strength  with  188.00  lb-in  Applied  Moment 
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Matrix  Failure  Compression  Strength  Increased  to  600%  of  Tensile  Strength  with  200  lb-in  Applied  Moment 


Figure  22.  Matrix  Failure  Compression  Strength  Increased  to  600% 


of  Tensile  Strength  with  200  lb-in  Applied  Moment 
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Matrix  Failure  Compression  Strength  Increased  to  600%  of  Tensile  Strength  with  235  lb-in  Applied  Moment 


B.  SQUARE  PLATES  SUBJECTED  TO  TRANSVERSE  LOADING 


A  square  plate  was  used  to  examine  the  damage  progression  that  occurs  with 
different  boundary  and  loading  conditions.  The  plate  had  the  length  to  thickness  ratio  of  50. 
A  quarter  model  was  used  due  to  plate  symmetry.  Boundary  conditions  were  either 
clamped  or  simply  supported  on  all  plate  edges.  The  quarter  model  would  then  have 
adjacent  edges  clamped  or  simply  supported  and  the  other  two  adjacent  edges  symmetric 
as  shown  in  Figure  24.  The  applied  load  was  either  a  concentrated  load  at  the  point  or  a 
distributed  load  over  the  xy  plane.  Both  loads  were  applied  in  a  positive  z  direction  (into 
the  plane). 
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Figure  24.  One  Quarter  Plate  Model 


Figure  26  through  Figure  3 1  show  the  resultant  damage  for  each  combination  of 
boundary  and  loading  conditions.  Under  all  conditions,  the  plate  could  sustain  a  higher 
distributed  load  without  failure  than  a  plate  with  a  concentrated  load.  Complete  matrix 
failure  occurs  for  even  a  small  load  for  the  present  composite  materials  whose  properties 
are  the  same  as  those  used  in  the  previous  section.  There  is  no  contribution  to  composite 
strength  by  the  present  matrix  material. 
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A  plate  with  the  clamped  boundary  condition  begins  failure  along  the  boundary 
edge  and  around  the  center  of  the  plate.  Damage  then  proceeds  between  the  two  as  shown 
in  Figure  25  through  Figure  28.  Simply  supported  plates  show  damage  initiating  around  the 
center  of  the  plate  and  propagating  outward. 

The  condition  of  composite  failure  for  this  analysis  is  defined  as  when  the 
composite  is  unable  to  support  an  applied  load  and  maintain  structural  integrity. 
Concentrated  loads  cause  the  composite  to  fail  more  readily  in  the  clamped  condition  than 
in  the  simply  supported  condition  as  indicated  by  comparing  Figure  25  and  Figure  26  with 
Figure  29  and  Figure  30. 
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Figure  25.  Clamped  Plate  with  a  Concentrated  Load  at  the  center 


□ 

1  1 

1  1 

u 

1  1 

1  1 

1  1 

1  1 

1  1 

1  1 

□ 

oj 

1  1 

1  1 

1  1 

1  1 

1  1 

1  1 

1  1 

□ 

1  1 

□ 

1  1 

1  1 

1  1 

1  1 

1  1 

1  1 

iM 

□' 

tassJ 

m 

■ 

1  1 

1  1 

□ 

m 

F  ,1 

m 

Hi 

s 

1  1 

m 

l<X^ 

kkl 

IkXI 

kkl 

f«l 

!«} 

1  1 

1  1 

kk3 

M 

M 

lOCI 

KH 

K><1 

m 

1  1 

m 

m 

IkkI 

kvi 

m 

kki 

m 

m 

m 

kki 

M 

HXI 

1 . 1 

hX1 

m 

m 

1  i 

m 

m 

ij 

m 

133 

Total  Load: 


312.52  lb 
241.00  lb 
125.00  lb 
25.00  lb 


Symmelnc  Boundary  fiber  Direction 

◄ - ► 


Figure  26.  Clamped  Plate  with  a  Concentrated  Load  at  the  center 
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Figure  27.  Clamped  Plate  with  a  Uniformly  Distributed  Load 
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Figure  28.  Clamped  Plate  with  a  Uniformly  Distributed  Load 
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Figure  29.  Simply-Supported  Plate  with  a  Concentrated  Load  at  the  center 
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Figure  30.  Simply-Supported  Plate  with  a  Concentrated  Load  at  the  center 
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Figure  31.  Simply-Supported  Plate  with  a  Uniformly  Distributed  Load 
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VI.  CONCLUSIONS  AND  RECOMMENDATION 


The  micro-macromechanical  analysis  allowed  for  the  study  of  damage  progression 
in  laminated  fibrous  composites  subjected  to  bending  loads.  By  first  applying  the 
micromechanical  model  developed  by  Kwon,  the  composite  properties  were  determined  by 
the  material  properties  of  the  fiber  and  matrix.  The  stresses  were  determined  at  the  fiber  and 
matrix  level  by  using  the  micro-macromechanical  approach.  Failure  analysis  then  could  be 
done  based  on  the  stresses  in  the  fiber  and  matrix.  The  strength  and  failure  mode  of  the 
composite  were  found  from  this  study.  Degradation  of  failed  fiber  and  matrix  material 
properties  was  included  in  the  program  such  that  damage  propagation  could  transfer  the 
load  from  a  load  from  a  damaged  section  to  neighboring  intact  sections  resulting  increased 
stresses  in  the  intact  area.  Applied  loads  were  increased  incrementally  in  the  model  until 
complete  composite  failure  occurred. 

Good  agreement  was  found  between  the  computer  model  and  experimental  results 
for  a  carbon/epoxy  laminated  composite  with  a  circular  hole  at  its  center.  The  parametric 
study  involving  the  plate  clearly  showed  the  dependance  of  the  composite  ultimate  strength 
on  the  fiber  strength.  As  the  fiber  strength  was  changed,  the  composite  strength  changed  in 
nearly  a  directly  proportional  amount. 

When  the  computer  model  was  used  to  study  a  square  plate  with  transverse  loading, 
the  failure  pattern  was  as  expected  for  such  conditions.  The  square  plate  study  also  showed 
that  a  greater  load  could  be  maintained  on  the  plate  with  less  damage  when  a  uniformly 
distributed  load  was  applied  instead  of  a  concentrated  load. 

It  is  suggested  that  this  study  be  extend  to  cover  cases  involving  multiple  holes  to 
determine  how  their  proximity  to  each  other  in  a  plate  influences  damage  progression. 
More  experimental  results  need  to  be  obtained  to  further  verify  the  damage  model. 
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APPENDIX:  FINITE  ELEMENT  MATRICES 
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